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We show that the superconducting 2SC phase at high density and normal chiraly broken quark 
phase at low density is separated by the mixed non-uniform phase along the baryon density line. 



I. INTRODUCTION 



In QCD with two flavors at high baryon density the superconducting phase is the lowest energy state [1,2]. In 
this phase u,d quarks with two different colors (say "red", "green") create Cooper pairs whereas the "blue" quarks 
remain free. The chiral symmetry 5(7(2) x SU(2) is restored and only gauge color group is broken. On the other 
hand at lower density one expects the state with broken chiral symmetry. This can be a hadronic phase but the free 
quark plasma is not ruled out. In such a situation the sequence of the phase transitions along the density line (at 
zero or very low temperature) follows the chain: hadronic phase/quark phase/superconducting phase. All the phase 
transitions are expected to be first order. In this letter we present an explicit model describing the transition between 
£T} ' the normal quark and the superconducting phases. It occurs that the transition actually goes through the mixed state 
t-H , which consists of the non-uniform chiral phase and superconducting phase. 

The non-uniform chiral phase had been introduced in the context of the non-relativistic nuclear physics [3] and then 
was extended for the systems of high density plasma of quarks [4,5]. In this new phase the chiral fields condense in 

■ both scalar and pseudoscalar channels creating a static chiral wave along one direction in space with the wavelength 
2iv/\q\ where cf is a wave vector. It was shown that in some temperature - density range this configuration has lower 
energy then the uniform quark phase. On the other hand we have learned from the studies [1,2] that at higher densities 
one should compare the non-uniform quark phase with the superconducting phase rather then with the uniform quark 

t-H , phase. As a result one can show that the phase diagram of strongly interacting matter is more complicated and 
C^l ' contains a region of the mixed state. 

The high density matter exists in the Neutron Stars. If the density is high enough then the superconducting phase 

■ resides in their cores. In this situation one can also expect the mixed phase to be present there - the mixed state 
CIh! which separates uniform chiraly broken phase from the superconducting phase. We cannot prove that this mixed 

state is the same as the one considered in this paper but we provide an explicit realization of such a scenario. 
' In the next section we introduce a thermodynamic potential describing all the interesting phases. In the third 
f~| . section we present the numerical analysis and finaly we discuss the results in the Conclusions. 



X 



II. MODEL OF CHIRAL AND 2SC PHASES 

Let us consider the NJL model [6]: 



H 



d 3 x ■ V - M7o )V> - G [OM) 2 + (^75tV>) 2 ] - G'(^T 2 t A C l5 ^(i>T 2 t A C 75 i>)} (1) 



where ip is the quark field, \i the quark chemical potential and the color, flavor and spinor indices are suppressed. 
The vector f is the isospin vector of Pauli matrices and t A are three color antisymmetric group generators. The 
coupling constant G describes the interaction in the isospin singlet, Lorentz scalar and the isospin triplet, Lorentz 
pseudoscalar, quark-antiquark channels whereas G' describes the interaction in the color 3, flavor singlet, Lorentz 
scalar diquark channel. Both couplings are related through the Fierz transformation but we treat them as the 
independent parameters. There is an additional parameter (or function) which regularize infinities in the model. In 
our case we use the simplest momentum cut-off A. Our main results are independent of this choice. We analize the 
four-fermion point interactions in the mean field approximation using the general anstaz: 

M 

(rjjiP) = -—cosq- x (2) 
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where the directions of symmetry breaking in the isospin and color spaces are chosen arbitrarily. The wave vector q 
describes the chiral wave along some arbitrary direction in space. In the limit of vanishing q we recover usual uniform 
chiral phase. The gap parameter A describes the superconducting phase. Let us notice that (2) is not the analog 
of the LOFF phase in the superconductivity [7,8] which arises due to the difference in the chemical potentials of the 
pairing species. In our case it is rather the dynamical interaction between the spin and isospin degrees of freedom 
which effectively lower the energy [4,5]. 

Using formulae (1), (2) in the mean-field approximation we get: 

/( - A A* M 2 IAI 2 1 

d 3 x -S7 + MU- M7 o)?/> - -^{tnt'C^ - — (V^C^V) + — + Llj-j (3) 

where the Dirac operator is now space dependent U — cos q-x + ^75X3 sin q ■ x. One can eliminate this dependence by 
the unitary transformation to the new fermi fields ip — U 1 ^ 2 ^ which gives the free energy in the form: 

/( 1 A A* M 2 IAI 2 1 

d 3 x W(-iS • V - -E ■ qr 3 + M 7o - n)$ - -(V^CW) 1 - — (V^CW) + ^g+^t\ ( 4 ) 

where S is the vector of Dirac spin matrices. The hamiltonian (4) describes two kinds of interactions: the "chiral" part 
and the "superconducting part". The chiral part mixes particles and antiparticles leaving color and isospin structure 
untouched. On the other hand the "superconducting" part mixes colors and isospins of particles. To be more precise, 
one can introduce the plane wave basis 

^ f d 3 k (->.-> -> - i 

= z2 / / 1 u s( k ) a i,s( k ) exp(-ikx) + v s (k)b> a ] s (k) exp(to) \ , (5) 

s=i,2 J {2ir)^2E(k) 

where E(k) = \Jk 2j r M 2 , u s , v s are Dirac bispinors, a l a s (k) (b l a s (k)) is an annihilation operator of quark (antiquark) 

of color a, flavor j, spin s, and momentum k, satisfying usual anticommutation relations. In this basis hamiltonian 
(4) creates the 24 x 24 matrix which decays into three independent diagonal 8x8 blocks. Two of them are related 
to pairing between "red", "green" u, d quarks and antiquarks whereas the third one describes the pairing between 
"blue" quarks and antiquarks (which do not create Cooper pairs). The hamiltonian (4) can be diagonalized with the 
positive energy eigenvalues : 

Ai,±(fc) = ^{E ± {k)-^) 2 + \A\ 2 , \ 2t± {k) = ^{E ± (k)+^) 2 + \A\ 2 (6) 

E±(k) = \jk 2 + M 2 + ^ ±yj 'q 2 M 2 + k ■ q 

for u, d, "red", "green" quarks and A3,±(fc) = E±(k) for "blue" u, d quarks. The negative eigenvalues has opposite 
signs. After diagonalization the hamiltonian (4) takes the general form Hmf = H + Hq. The first term is a diagonal 
operator and the second term defines the energy of the ground state as a function of chemical potential. The vacuum 
expectation value (Hmf) = H and 

IAI 2 M 2 „^ f A d 3 k /x , nV ^ f A d 3 k , f d 3 k 

The first integral describes the contribution from the quarks which create the chiral condensate as well as Cooper 
pairs. The factor 2 is the number of colors minus one ("red", "green"). The parameter A reminds of the momentum 
cut-off of the divergent integral. The last two terms are connected to "blue" quarks which build only chiral condensate. 
The factor of 2, this time, describes the number of flavors. Let us rewrite (7) in the form: 

|A| 2 M 2 s—^ f A d 3 k n ^ \ f d 3 k „ nl , f d 3 k „ .1 /n , 
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which distinguishes explicitly between the vacuum contribution - the first integral, and the finite density contributions 
- the last two integrals. The vacuum contribution can be expanded in the power of the wave vector q (nothing more 
but the gradient expansion in mesonic fields) and we arrive at the final result: 



A| 2 M 2 M 2 F 2 q 2 _ f A d 3 k 
!_ j I " —19 / 

±G' AG 2M 2 J 



r A H 3 k 



The coefficient at q 2 is model independent and it is related to the pion decay constant with its vacuum value = 93 
MeV. M is the constituent quark mass at zero density. Let us notice that in the limit of vanishing |A| we recover the 
chiral hamiltonian [5] whereas in the limit of vanishing q we can find the hamiltonian describing Cooper pairing [6] . 
Let us also mention that from the formula (9) one can easily guess the linear sigma model version of the hamiltonian. 
Indeed the NJL and linear sigma models at the mean-field level are different only with respect to the description of 
the vacuum contribution to the energy. Thus instead of the formula (9) one can write in the linear sigma model: 

h* |A|2 + M2 + M2p2 ^ + m " ( m2 fA 2 rim 



where g is dimensionless coupling constant and m a is a mass of sigma meson. However in the numerical calculation 
we only use the formula (9). 



III. CHIRAL/2SC PHASE TRANSITION 

The model parameters can be fitted to the values of the chiral condensate and the pion decay constant which give 
G = 5.01 GeV~ 2 and the cut-off A = 0.65 GeV [9]. For this set of parameters the chiral phase transition at q = 
and |A| = takes place at the chemical potential fi — 0.316 GeV and the value of the constituent quark mass at zero 
density is M = 0.301 MeV. The coupling constant G' we treat as the additional model parameter. The main results 
which are robust against the change in G' are: 

• there is the first order phase transition from the homogeneous chiraly broken phase M = Mo, q = 0, |A| = to 
the mixed phase M ^ M ^ 0, q ^ 0, |A| ^ 0. 

• by increasing chemical potential one decreases constituent mass M and increases q and |A|. 

• there is additional first order phase transition from the mixed phase to the 2SC phase (M = 0, q = 0, |A| ^ 0). 

From the above points one can see that the 2SC phase is separated from the chiraly broken phase by the the non- 
uniform mixed state. This conclusion is still valid for non-zero temperatures. We discuss this feature in the last section. 
On the other hand the numerical details depend on the choice of G'. In particular for the range G/2 < G' < G we 
can find: 

• The transition to the mixed phase appears in the band of 10 MeV around [i = 0.29 MeV (smaller G' larger 
critical chemical potential) . 

• The values of the gap parameter |A| can vary by a factor of two (smaller G' larger |A|). 

• The values of the wave vector |<f| can change by 40-50 per cent (smaller G' larger \q\). 

• The range of the mixed phase along /x-axis can change from 10 to 100 MeV (smaller G' longer the range). 

In the Fig. 1 we present the values of M, |A| for G' = G/2 as a function of chemical potential /i. The first 
order phase transition takes place at fi = 0.301 GeV. The constituent mass drops to around half of its vacuum value 
(Fig. la), the wave vector |<f| = 0.34 GeV (Fig. lb) and |A| = 0.004 GeV (Fig. lc). The mixed phase exists up to 
the second critical point /i = 0.396 GeV where another first order phase transition appears. The values of M, q drop 
to zero and the value of the gap A| = 0.039 GeV. For larger fj, we enter the homogeneous 2SC phase. 
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FIG. 1. Dependence of M, \q\ and A in GeV on the quark chemical potential n in GeV. 



IV. CONCLUSIONS 



In this letter we showed that the high density superconducting phase 2SC is separated from the low density chiraly 
broken phase by the region of the mixed non-homogeneous phase. This mixed state consists of the chiral waves with 
non-zero expectation values in the scalar and pseudoscalar channels (tpip), (ipij^Tatp) described by the wave vector q 
which points at arbitrary direction is space. In addition there is a non-zero value of the superconducting gap parameter 
A|. We cannot prove that this mixed state is the lowest energy state however we showed that its energy is lower then 
that of the homogeneous phases. Similar questions was already addressed in the context of the Overhauser effect [10] 
versus BCS pairing. It was shown that both types of pairing can be competitive [11,12] with each other. However 
it was also shown that within the NJL-type of models BCS state is favoured over the Overhauser state [12] which 
suggests that the non-homogeneous phase (2) is lower then BCS and Overhauser phases in the considered density 
range. Nevertheless one has to keep in mind that other approaches (instanton models) favours Overhauser effect over 
BCS pairing [11,12] thus one can not drew the final conclusion and more work is required. We can then conclude that 
the superconducting phase and the normal phase are separated by some kind of crystal state. Let us stress that the 
crystal structure (2) is dynamical in origin, not " kincmatical" in a sense of the LOFF state [7,8]. 

Our results are robust against the reasonable changes in the model parameters. What is more the same results one 
can derive not only in the NJL model but also in the linear sigma model. Thus also in that sense the results are model 
independent. The numerical details, like the location of the critical points, the values of the constituent mass, wave 
vector and superconducting gap, depends on the choice of the parameters. Although we perform our calculations for 
massless quarks the introduction of the small current quark mass does not change the main conclusions. It merely 
distorts the shape of the chiral wave as was shown in [13] for one dimensional system. One can also expect that the 
mixed phase exists for some range of temperatures (below T ~ 50 MeV where superconductor melts) and densities. 
However the detailed picture of the phase diagram remain to be done. 

The real QCD consists of two light flavors and one heavier, strange quark. In this situation one has to consider 
the three flavors NJL model. However we expect that our main points remain the same. If one considers the 
superconducting phase inside the Neutron Stars then the question of the color and charge neutrality rises [14]. It 
was shown by the explicit calculation [15] that at non-zero lepton chemical potential the 2SC phase borders with the 
normal quark phase at lower densities and only at higher densities it changes into the CFL phase. Thus even in the 
case of the three flavors QCD we can still have quark matter/2SC boundary in the Neutron Stars. This boundary 
is not just a line in the phase diagram but the mixed state region as we have shown. What kind of the crystal or 
non-uniformity is realized in the three flavors QCD has not been discoverd so far. 
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